Systematic analytical approach to correlation functions of resonances in quantum 

chaotic scattering 
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We solve the problem of resonance statistics in systems with broken time-reversal invariance by 
deriving the joint probability density of all resonances in the framework of a random matrix approach 
and calculating explicitly all n-point correlation functions in the complex plane. As a by-product, 
we establish the Ginibre-like statistics of resonances for many open channels. Our method is a 
combination of Itzykson-Zuber integration and a variant of nonlinear a— model and can be applied 
when the use of orthogonal polynomials is problematic. 



PACS numbers: 05.45.+b 

As is well-known, universal statistical properties of 
bound states in the regime of quantum chaos can be de- 
scribed in the framework of the random matrix approach 
The relevant methods adjusting random matrix de- 
scription to the case of resonance scattering in open quan- 
tum systems are very well-known since the pioneering 
work by the Heidelberg group j3), see the review jjj for 
a thorough discussion of recent developments. 

One of the basic concepts in chaotic quantum scatter- 
ing is the notion of resonances. Resonances are long-lived 
intermediate states to which bound states of a "closed" 
system are converted due to coupling to continua. On the 
formal level the resonances show up as poles of the MxM 
scattering matrix S a b{E). The dimension of this matrix, 
M, equals the number of open channels in a given inter- 
val of energies. The poles of S a b{E) occur at complex 
energies £u = E k — |r fe , where E k is called the position 
and Tfc the width of the corresponding resonance. Recent 
advances in computational techniques made available res- 
onance patterns with high accuracy for realistic models 
of atomic and molecular chaotic systems ||, as well as 
for quantum billiards and other models related to chaotic 
scattering ||. 

In the framework of the random matrix approach the 
S-matrix poles (resonances) are just the complex eigen- 
values of an effective random matrix Hamiltonian H e s — 
H — iT. Here H is a random self-adjoint matrix of a 
large dimension N describing the statistical properties 
of the closed counterpart of the scattering system un- 
der consideration. Depending on presence or absence of 
the time-reversal invariance H has to be chosen as a real 
symmetric or complex Hermitian one, respectively jUI). 
The N x N matrix V serves the purpose of describing 
transitions from the states described by H to the outer 
world via M open channels. It is related to the N x M 
matrix W of transition amplitudes in the following way: 
T = nWW^ . Such a form of T is actually dictated by the 
requirement of the S— matrix unitarity and ensures that 
all 5— matrix poles lie in the lower half-plane of complex 



energies, as required by causality. It is evident that the 
rank of T is M. In practice, the most interesting case is 
that of few open channels: N M ~ 1. In this case 
the width 1^ of a typical resonance is comparable with 
the mean separation A between neighboring resonances 
along the real axis and statistical properties of resonances 
are expected to be universal . 

Despite quite substantial efforts (?],||,|],|llj our actual 
knowledge of S-matrix poles statistics for few-channel 
scattering is still quite restricted. Among the facts estab- 
lished analytically beyond perturbation theory one can 
mention (i) the density of joint distribution of all res- 
onances for the system with a single open channel and 
gaussian-distributed transition amplitudes W and (ii) 
the mean density of S-matrix poles for arbitrary M<iV 
§§, as well as for M~N @. 

At the same time, the most interesting and difficult 
question of correlations between resonances in the com- 
plex plane resisted systematic analytical investigations. 
As an attempt to get an insight into the problem a non- 
trivial integral relation satisfied by the lowest (two-point) 
correlation function of complex eigenvalues for T > was 
derived recently in Q . Using that relation it turned out 
to be possible to put forward a conjecture on the an- 
alytic structure of the correlation functions for systems 
with broken time-reversal invariance. Unfortunately, the 
above mentioned relation neither fixed the lowest cor- 
relation function in a unique way nor provided a direct 
information on higher correlation functions. 

The goal of the present paper is to develop a regular 
analytical approach to the statistical properties of reso- 
nances for systems with broken time-reversal invariance. 
To this end, we first derive the joint probability density of 
all resonances for arbitrary number of open channels M. 
Then we reduce the problem of extracting the n— point 
correlation functions in the limit n, M-fixed, N — > oo to 
averaging a certain product of 2n determinants. Finally, 
the latter is evaluated with a method combining a map- 
ping to a fermionic version of a nonlinear ct— model with 
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the Itzykson-Zuber integration. As a result, we prove the 
validity of the conjecture put forward in JTTQ - 

We consider an ensemble of random N x N complex 
matrices J = H + iT |jl3), where H is N x N matrix 
taken from a Gaussian Unitary Ensemble (GUE) of Her- 
mitian matrices with the probability density V(H) oc 
cxp (— y tr H 2 ) and T is a fixed nonnegative one: T > 0. 
The probability density function in our ensemble can be 
written in the form 



V(J) oc exp 



N fJ + ,P 

— tr 

2 I 2 



s r- 



j-jt 

2i 



(1) 



Here and henceforth we do not specify the multiplica- 
tive constants when dealing with probability densities 
and correlation functions since they can be always found 
from the normalization condition. 

Eq. (|l|) can be used to obtain the density of joint dis- 
tribution of eigenvalues by integrating V( J) over the de- 
grees of freedom that are complementary to the eigen- 
values of J. This can be done following Dyson's method 
(see [^|,^0|). Neglecting matrices with repeated eigenval- 
ues, one transforms J to triangular form: J = U{Z + 
R)U\ where U is a unitary matrix, R is strictly upper- 
triangular and Z = diagjzi, . . . zn} is the diagonal ma- 
trix of complex eigenvalues. The Jacobian of the trans- 
formation from J to {Z,U,R) is |A(Z)| 2 where A(Z) = 
T\i<j<k<N^ z i ~ Zk ) IS ^he Vandermonde determinant. 
To perform the integration over R it is technically con- 
venient to use a Fourier-integral representation for the 
5— function in Eq. (Q). This reduces the corresponding 
integral to a Gaussian one and after algebraic manipula- 
tions the resulting expression is 



P N ,r(Z)cxe-* RctrZ2 - 



■ trT* 



|A(Z)| 2 Q(ImZ), (2) 



where Q(ImZ) is the remaining integral 

N 

Q(Im Z) = J [dU] J] 5 (Im z x - {U^TU) U ) 



(3) 



over the unitary group U(N), [dU] being the Haar mea- 
sure. Due to the specific structure of the matrices J their 
eigenvalues lie in the upper part of the complex plane and 
in all formulae below Im Zj > for all j JT3J . 

To proceed further we need to integrate over U. Again 
it is convenient to use the Fourier-integral representation 
for the S— functions in Eq. (||): 



Q(ImZ) 



dK 



J, Im tr KZ 



(27T) 



N 



[dU] e ~ itlKU ' ru , (4) 



where the first integration is over all real diagonal matri- 
ces K of dimension TV, dK being dk\ . . . dkjq . 

When the eigenvalues of T are all distinct the integra- 
tion over U can be performed using the famous Itzykson- 
Zuber-Harish-Chandra (IZHC) formula [Q. We, how- 
ever, are mostly interested in the case when T has a small 



rank M <C N, i.e. it has only M nonzero eigenvalues 
which we denote by 71 , ... , jm ■ This limit of highly de- 
generate eigenvalues is difficult to perform in the original 
IZHC formula. Nevertheless the difficulty can be circum- 
vented and the result is as follows: 



QQmZ) 



det 



M-N 



7 



A(7) jr 

N M 



dA det[ e - l7iA " 



(5) 



x A ( A )nE 

diag(Ai 



iX lrl IT: 



, Am) and 7 = diag(7i, 



,7m)- 



where A = 

Eqs. (||) and ([5|) give an explicit representation for the 
joint probability density of N resonances Zi in the com- 
plex plane. As such, they constitute one of the main 
results of the present paper and provide the basis for cal- 
culating the n-eigenvalue correlation functions 

JV! f 

^( z )= (N-n)\ J dwP ^r(z,w), (6) 

where, for the sake of brevity, we decompose Z = 
diag(z, w) with z= (zi, ...,z n ) and w = (w b ... , w N - n ), 
dw = Y[J = i n dRewj dlmwj, identifying Wk = z n+ k- 

In what follows we will calculate i?„(z) for arbitrary 
fixed n and M in the limit N — > 00. On the first stage 
we will replace the integration over w in (^|) by averag- 
ing over the ensemble of non-Hermitian random matrices 
Jn-u(j) = Hjf-n + *r, with -Hjv-n being a GUE ma- 
trix of the reduced size (N — n) x (N — n). This step 
involves cumbersome algebraic manipulations and will 
be presented in full details elsewhere. Here we outline 
the ideas on the simplest, still nontrivial example of one- 
channel systems (M = 1 ). In this case T has only one 
nonzero eigenvalue 7 and the A— integration in Eq. (|^) 
can be explicitly performed yielding 

|A(Z)| 2 
Pjv, 7 (Z)cx _ w=r 

1 " 3=1 



i V 

1(7 -^Im; 



(7) 



Y] Im z j , we write 



Introducing the notation 7 = 7 

the 5-function in Eq. (Q) as S (^y — Im w/J . Now 

we replace N in the exponent of Eq. ([7]) by N — n (this 
act is justified by the limit N — > 00). With the resulting 
relation in hand one readily obtains that 

C7O) 



R n (z) cx 



|A(z)| 



2 e- ! 



N—n—l 

e~ 



-(7 2 -7 2 ) 



(8) 



where 



N—n n 



c^z) = /dw p w _„,^(w) n n 



1=1 ] = 1 



Yl det [Zj - JN-n(l)} 



Wl\ 



GUE 



(9) 
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In the limit when N — > oo and M is finite (in partic- 
ular, for the present case M = 1), the imaginary part 
of almost all eigenvalues of J is of the order -k <C 7 m 
jl4j], hence so is 7 — 7. Therefore one can reinstate 7 in 
place of 7 in the determinants in Eq. (0) . On the other 
hand, rescaling the imaginary parts yj — NlvaZj, one 
finds that 



N-n-l 



in the limit n <C iV — > 00, with <? = ^(7 + 7 _1 ). 

Essentially similar manipulations can be performed for 
arbitrary fixed number of open channels M. In the limit 
n, M <C N — > 00 we arrive at the following general rep- 
resentation of the correlation functions: 



i?„(z) cx 



det" 7 

n M 



a( z ; 



I 2 



(10) 



j=l m 



1 ris^m(5m 9s) 



where g m = 5(7™ + jj-). 

Thus, the problem amounts to evaluation of the corre- 
lation function of the determinants in Eq. (||) Q . To pro- 
ceed, we first write each of the determinants as a Gaus- 
sian integral over a set of Grassmann variables. When 
this is done, the GUE average becomes trivial and yields 
terms quartic with respect to the Grassmannians. These 
terms can be further traded for an auxilliary integration 
over a Hermitian matrix S of the size 2n x 2n (the so- 
called Hubbard-Stratonovich transformation). Then the 
integration over the Grassmann fields is trivially per- 
formed and yields again a determinant. As the result, 
we arrive at the following expression: 



(11) 



M 



J\ det [l 2n + «7m L 2n(Z 2 „ - iS) *] 



where we have introduced the diagonal matrices Z 2n = 
diag(z,z t ) and L 2 „ = diag(l n , -1„). 

Let us now recall that nontrivial eigenvalue correla- 
tions are expected to occur on the scale when the 
eigenvalues are separated by distances comparable with 
the mean eigenvalue separation for GUE matrices H, 



the latter being of the order (N 



with our choice 



of V(H). Accordingly, it is convenient to separate the 
"center of mass" coordinate x = ^Y^j=i^ ez j so that 
Zj = x + ffl n j where both the real and imaginary parts 
of Sj are of the order of 1 in the limit when N — > 00 and 
M is fixed. In this limit R n (z) is effectively a function of 
z (x is fixed) which we are going to calculate. 

To evaluate the integral in ( pi] ) let us first diago- 
nalize S: S = U2nSU 2 ~ n , where J7 2n € U(2n) and 



S = diag(eri, er 2n ). Then, keeping only the leading 
terms in the limit N — > 00, we obtain: 

O r {z) = /d£A 2 (£)e (C(S)) u(2n) 



(12) 



where 

(C(S)) u{2n) = J [dU 2n ] e 

M 



-tt[Z 2n {xX 2n -iSy 1 ] 



Y[ det [l 2n + ij m h 2n (x±2n - iS) x ] . 



The form of the integrand in ( |12[ ) suggests exploit- 
ing the saddle-point method in the integral over er^, 
k = 1, ...,2n. Altogether there are 2 2n saddle-points: 
cr^' = — |(a; + iekV^ — x 2 ), where = ±1. The leading 
order contribution comes from integration around those 
saddle-points where exactly n parameters equal 1 (the 
rest being equal -1). All other choices can be neglected 
as they lead to lower order terms. This is because of 
the presence of the Vandermonde determinant in the in- 
tegrand. At the same time, all relevant saddle-points 
produce the same contribution and we obtain that: 

C T (z) cx e^^U^^c^l) (13) 

where z = [zi, . . . , z n ), Zj = N(zj — x), j = 1, . . . ,n, and 



C;(i) = J [CI0 2 „J r 
M 



— i7Tl/(x) tr Z2tiQ2t7 



(14) 



JJdct 



llcX, 



i2n H ^— L 2n + 7rzy(a;)7 c L 2 „Q 2 „ 



In (|1J) Q 2n = U 2n L 2n /7 2 „ , the integration is over the 
coset space U(2n)/U(n)®U(n), and the symbol v(x) 
stands for the semicircular density of real eigenvalues of 
the matrices H, v(x) — ^:V4 — x 2 . 

Thus, the problem reduces to evaluation of an integral 
over a coset space. This type of integrals is known in the 
literature under the name of zero-dimensional nonlinear 
a— models and our consideration enjoy many useful in- 
sights, see e.g. fig ]. In particular, the following polar 
parametrization proves to be the most effective: 



U 2 „ = 




cos 



i<t> 



sin^ 



sin ip — cos ip 




where Ua,r & U(n), if) = diag(^i, tp n ), <j> = 
diag(0i, (f> n ) and < 4>k,ipk < 2ir. The corresponding 
measure [dQ 2 „] is proportional to 



[ci[/,4] [dL7?J Y[ d(j>jdipj sin 2^- JJ | cos 2^— cos 2ipk 



l<l<k<n 
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The use of such a parametrization makes the integra- 
tion especially simple because of the determinant factor 
being independent of the unitary matrices Ua,r- As a re- 
sult, these matrices appear only in the exponential factor 
and the corresponding integrals can be evaluated accord- 
ing to the IZHC formula |l2|]. Passing to the variables 
Xk = cos 2ipk we find: 



det" 7 

WW 

x det 
n! det™ 7 



Y\_ d\j jj Gm(Aj 



,;=i 



3=1 



where 



1 



,1/ 



Gat (A) = [I I.9m + 7Ti/(a;)A] 



m— 1 



Combining this with Eqs. ( jToj ) and (|i~3| ) and restoring the 
normalization we finally see that the correlation func- 
tions, in the limit N — > 00, have the following simple 
structure: 



N 2r, 



-R„ [x + —,...,x+ =&et[K(z ,z* k )] n jk=l 



N 



where the kernel K(Zj, z k ) is given by 
K{z 1 ,~4) = F 1/2 {z l )F 1 ' 2 {z* 2 ) 



(15) 



with F(~z) = E^ =1 



n, 



This is equivalent to 



the form conjectured in |11|. 

One of the important physical limits of the scatter- 
ing system is the case of many equivalent open channels: 
<Jm = g for all m and M 3> g. Resonances in that case 
form a dense cloud in the complex plane characterized by 
a mean density p(z) inside the cloud. This fact and ex- 
pression for p(z) were found in Using our formulas 
derived above we are able to show that the statistics of 
the resonances in that case is determined by a Ginibre- 
like kernel: 

\K(zi,z 2 )\ = p(z)enp--np(z)\zi - z 2 \ 2 (16) 

with z = (zi + z 2 )/2, generalizing a classical result by 
Ginibre [[TtJ to the case of a nonuniform (i.e. position- 
dependent) mean density of complex eigenvalues p(z). 
As such, it has a good chance to be universally valid for 
strongly non-Hermitian random matrices. 

In conclusion, we considered a non-Hermitian random 
matrix model of chaotic quantum scattering. We started 



with deriving the joint probability density of all com- 
plex eigenvalues describing S-matrix poles ( resonances) 
in chaotic systems with broken time-reversal invariance. 
Then we found a way to extract all correlation functions 
of the resonances in complex plane. As a by-product, we 
established the Ginibre-like statistics of resonances for 
many open channels. 
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